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Motivation 



How well is the vehicle’s aerodynamic performance estimated? 

Is the mesh appropriate for every flow condition and vehicle configuration? 



Objectives 

Toward automation of CFD analysis 



Handle complex geometry problems 

Control discretization errors via solution-adaptive 

mesh refinement 

Focus on aerodynamic databases of parametric 
and optimization studies 

1 . Accuracy: satisfy prescribed error bounds 

2. Robustness and speed: may require over 10 5 
mesh generations 

3. Automation: avoid user supervision 

Obtain “expert meshes” independent of user skill 
Run every case adaptively in production settings 


Approach 


1 . Embedded-boundary Cartesian mesh method (1 990’s) 

• Arbitrarily complex domains, efficient and accurate 

• Irregularity confined to body intersecting cells 




2. Incremental strategy for h-refinement of 
nested Cartesian meshes (2002) 

• Fast local re-meshing of flagged cells 

• Guaranteed reliability 

• Early work used feature detection and t- 
extrapolation 


ort3D 


3. Adjoint-weighted residual error estimates (2007) 

• Mesh enrichment targets output functionals 

• Functional error-bound estimates 

• Implementation exploits nesting of Cartesian 
meshes for fast interpolation 





Functional (Cl, Cd, . 


Numerical Error 
Uniform Mesh Refinement 




Number of Cells 


> Exact Functional: J 

• Numerical solution on a mesh 
with cell-size H gives 
approximate functional: 

J(Ujf) 

• Goal is to estimate functional 
error: 

J(U) = J(Uff) + E 

• Express the error as a function 
of the flow solution 

E = f( U fl ) 


Functional (Cl, Cd, . 



Discrete Estimate of Numerical Error 



Number of Cells 


Consider a simpler problem of 
computing relative error: 


J( Ufc) = J(Vh) + e 


For second-order accurate 
spatial discretization and cell- 
size in the asymptotic range, the 
functional error is: 

1 1 

E = e + -e - b —^e J \ 

4 4 2 



We will use an adjoint solution 
on mesh H to estimate 


e = f (Ufl.fe) 


Adjoint Error Estimates 


• Consider a functional J (Ujj) computed from the solution of Euler 
equations discretized on an affordable mesh with cell-size H: 


R(Vh) = 0 



• In addition, consider an embedded mesh with cell-size h obtained via 
uniform refinement of the baseline mesh 

• We seek to compute the error relative to the embedded mesh without 
solving the problem on the fine mesh 

e = | J(U„) - J(U")| 


Venditti & Darmofal, 2002 




Estimate of functional on embedded mesh is obtained from Taylor 
series expansions about the coarse mesh solution 


J{ Ufc) 
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These equations are combined to give 


J(U ft )« J(Uf ) - ^ R(U£ ) 

where 'p satisfies the adjoint equation 
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• Estimate of functional on embedded mesh is obtained from Taylor 
series expansions about the coarse mesh solution 

J(U h ) « J( u") + - Uf ) 


R(U h) 


0~R(U") + 


^R(uf) 

<9u h 




These equations are combined to give 


J(U ft ) * J( Uf ) - 


where 'p satisfies the adjoint equation 






Adjoints provide 
a weighting on 
residual errors 


Adjoint Correction and Error Bound 



• Since the adjoint solution is not known on the embedded mesh, 
we use an approximate solution from the coarse mesh to obtain 

j{ u h ) « A u%) - {< ) T R(Uf ) - - ^) T R(Uf ) 

Adjoint Correction Remaining Error 

rx zj 

• XJ h ,ip h denote reconstructed solutions lifted from coarse mesh 
to embedded mesh. We use linear interpolation 

• 'iph is unknown. We approximate it with a 
quadratic interpolant 




Adjoint Correction 
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Number of Cells 


Predict functional on a fine 
mesh with cell-size h from a 
coarse mesh solution with cell 
size H 


U H 



Uh 








Adjoint Correction 



(V>/i 


<) T R(Uf) 



Predict functional on a fine 
mesh with cell-size h from a 
coarse mesh solution with cell 
size H 


U H 



Uh 








Error Bound Estimate 


J(U„)«J(U?) 


«>f) T R(Uf) 



• Bound on remaining error in 
each coarse cell k 

efc = El^Q-^) T «(U L ) 

N 

• Net functional error E — E 

k= 0 
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• Given a user specified tolerance TOL, termination criterion is 
satisfied when E < TOL 




Refinement Parameter 


• Define maximum allowable error level in each 
coarse cell via equidistribution: t = TOL / N 

• Refinement parameter in each cell is given by r k 

• Refine cells for which r & > A 
where X > 1 is a global threshold factor 


Error Histograms 




t 
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Perfect 

Mesh 

► 


Cell-wise Error 





Results 

Focus on Applications 



Part A. Accuracy 

• Launch Abort Vehicle with jets - uniform mesh refinement study 

Part B. Efficiency 

• Sonic-boom signature test case - computational cost summary 

Part C. Databases 

• Nozzle-Guide-Vane Missile 

• Launch Abort Vehicle with Jettison Motor plumes 


Launch Abort Vehicle 
with ACM Jets 



Mach -0.75 
6 sec 

-5300 ft altitude, -3300 ft downrange 


Mach -0.75 




Mach -0.3 
15 sec 


ACM trims LAV at -0 deg alpha during 
ACM coast phase. Slows down from 
Mach 0.8 to 0.2 in -10 seconds and 
drops in qbarfrom 800 to 100 psf 



Mach -0.2 
20 sec 


ACM holds LAV at -25 deg 
alpha during AM burnout. 



Mach -0.3 
1 sec 



ACM + Canards turn LAV to 
heat shield forward 


Pad-abort mission profile 


ACM initially pushes 
LAV to -25 deg alpha 
during AM burn 



8 Abort Control 4 Abort 
Motors (AC Ms) Motors (AMs) 


Mach -0.2 
21 sec 


LAS is jettisoned with ACM 
still burning. CM in free flight 
for -3 seconds until drogue 
chute deployment. 




Mach 0 
0 sec 


Ignition 



c 


4 Separation 
Motors (SMs) 


C.G. 



Launch Abort Vehicle 
with ACM Jets 



Mach -0.75 


Mach -0.75 
6 sec 

-5300 ft altitude, -3300 ft downrange 


Mach -0.75 




Mach -0.3 
sec 


A 15 s< 


Mach -0.2 
20 sec 


/// 


ACM trims LAV at -0 deg alpha during 
ACM coast phase. Slows down from 
Mach 0.8 to 0.2 in -10 seconds and 
drops in qbarfrom 800 to 100 psf 


X 


ACM holds LAV at -25 deg 
alpha during AM burnout. 


c 


Mach -0.3 
1 sec 


ACM + Canards turn LAV to 
heat shield forward 



ACM initially pushes 
LAV to -25 deg alpha 
during AM burn 


8 Abort Control 
Motors (ACMs) 


4 Abort 
Motors (AMs) 


Mach -0.2 
21 sec 


LAS is jettisoned with ACM 
still burning. CM in free flight 
for -3 seconds until drogue 
chute deployment. 


Mach 0 
0 sec 


Ignition 


'i 



4 Separation 
Motors (SMs) 



C.G. 
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Launch Abort Vehicle with ACM Jets 
Problem Setup 



Examine aerodynamic performance 
with ACM jets (AIAA 2008-1281) 

Selected case: Moo = 4, a = 20°, due to 
significant plume penetration 

Power boundary conditions applied at 
plenum face (assumes perfect gas) 

Thrust setting 

0 ibf 

0 Ibf 0 Ibf 


180 Ibf 


80 Ibf 



1819 Ibf 


Nozzle geometry 


2557 Ibf 



Launch Abort Vehicle with ACM Jets 
Functional 



• Functional: Cn+0.4Ca 

• Consider two cases: 

► Case A: Functional defined over 
nose-cone surface only, T0l_=0.0005 

- Accuracy verification: Uniform 
mesh refinement study 

- Take advantage of supersonic 
freestream conditions to limit 
refinement to nose-cone region 

► Case B: Functional defined over 
entire vehicle, T0l_=0.006 

- Typical engineering database case 



Launch Abort Vehicle 

Initial mesh and solution on symmetry plane 




Mach contours, M« = 4, a = 20° 



Case A: Finest Mesh of Uniform Refinement 
Study (Side-view, 75M Cells) 




Case A: Finest Mesh of Uniform Refinement 
Study (Side-view, 75M Cells) 



Case A: Finest Mesh of Uniform Refinement 
Study (Front-view, 75M Cells) 




Case A: Adapted Mesh (Side-view) 
1 4 Cycles; 31 Ok Cells 




Functional 


Case A: Functional and Error Convergence 




• Difference in functional values is below 0.05% on finest mesh 

• Two orders-of magnitude savings in total number of cells 

• Adaptive computation required just 9 minutes of wall-clock time on an 8-core Intel 
Xeon desktop 


Case B 

Final mesh and solution on symmetry plane 




15 adaptations, Mach contours, M® = 4, cx = 20' 



Case B 

Type IV Shock-Shock Interference 
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Close-up view of lower surface ACM 
|1 and plume, colored by Mach number 


Case B 

Final mesh at various x-stations 




15 adaptations, M® = 4, a = 20° 




Case B 

Front view of plumes on final mesh 



Launch Abort Vehicle 

Plume Visualization on Final Mesh (-7.7M cells) 



Moo = 4, a = 20' 



Side-view: plumes shown as iso-surfaces of total 
temperature colored by Mach number. Also shown are 
Mach number on symmetry plane and C p shading on body 

Main jet interaction occurs as lower-surface plume strikes 
sides of the boost protective cover. 

Largest errors in functional are near edges of main plume 
near the ACMs 


Launch Abort Vehicle 

Bottom view of plumes on final mesh 



Paths of three main plumes from the lower surface ACMs to the heat-shield. Main jet 
splits into two plumes that contact the aft region of the vehicle 


Launch Abort Vehicle 
Functional Convergence 


• Determination of plume paths and 
appropriate refinement of plume edges 
is not possible a-priori, yet these 
features determine the “aerodynamic 
shape” of the vehicle 


• Adjoint error analysis identifies regions 
where jet interaction effects are 
important for the computation of 
aerodynamic coefficients and triggers 
mesh refinement 


• Functional convergence settles down 
at ~1M cells, however, additional 
research is required to improve 
estimates of the error-bound 


C..+0.4C 




Number of Cells 


Results 

Focus on Applications 



Part A. Accuracy 

• Launch Abort Vehicle with jets - uniform mesh refinement study 
\ 

Part B. Efficiency 

• Sonic-boom signature test case - computational cost summary 
) 

Part C. Databases 

• Nozzle-Guide-Vane Missile 

• Launch Abort Vehicle with Jettison Motor plumes 


69° Swept Delta Wing-Body 


• NASA TN D-71 60 

► Moo = 1.68 

► a = 4.74° 

► Sensor offset, h/L = 3.6 & 

{0.2, 0.4, 0.8, 1.2, 2.0, 2.8} 

0.2 L 

• Initial mesh ~ 22 k cells 






69° Swept Delta Wing-Body 


• NASA TN D-71 60 

► Moo = 1.68 

► a = 4.74° 




69° Swept Delta Wing-Body 


NASA TN D-71 60 

► Moo = 1.68 

► a = 4.74 



Isobars 


69° Swept Delta Wing-Body 


NASA TN D-71 60 

► Moo = 1.68 

► a = 4.74° 



0.4 L 
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69° Swept Delta Wing-Body 


• NASA TN D-71 60 

► Moo = 1.68 

► a = 4.74° 





69° Swept Delta 


NASA TN D-71 60 

► Moo = 1.68 


► a = 4.74° 

► h/L ={.2, .4, .8, 1.2, 2.0, 2.8, 3.6} 

Simulation performed on 
desktop workstation 

► Dual quad-core (8 cores) 

► Intel Xeon, 3.2Ghz 

► 16Gb memory 



Adjoint Solv 
38 % 


Flow Solves 
49 % 


Mesh 
Adaptatio 
13 % 


Total simulation time 53 mins, 
(all adaptations & mesh gen) 



Total = 53 mins. 


Results 

Focus on Applications 



Part A. Accuracy 

• Launch Abort Vehicle with jets - uniform mesh refinement study 

Part B. Efficiency 

• Sonic-boom signature test case - computational cost summary 


Part C. Databases 

\ 

• Nozzle-Guide-Vane Missile 


• Launch Abort Vehicle with Jettison Motor plumes 



J 


Error Controlled Aero Database 


• Realistically complex model with plenum, guide vanes, etc. 

• Perform (data blind) aero analysis over range of operating conditions 

• Moo = {0.5, 0.7, 0.9, 1 .1 , 1 .3, 1 .6, 2.0} 

• a = {0°, 5°, 10°, 15°, 20°} = 35 cases total 

• Output functional: J = Cn + 0.1 Ca, tol = 0.05 


• Starting mesh has -8000 cells 




Error Controlled Aero Databas 


• Realistically complex model with plen 

• Perform (data blind) aero analysis ove 

• Moo = {0.5, 0.7, 0.9, 1 .1 , 1 .3, 1 .6, 2 .( 

• a = {0°, 5°, 10°, 15°, 20°} = 35 cases total 

• Output functional: J = Cn + 0.1 Ca, tol = 0.05 


• Starting mesh has -8000 cells 



Example Case 


M m = 0 . 9 , (X = 10 ° 



Example Case 

Mod = 0.9, a = 10° 




Error Controlled Aero Database 



(X = 0 ° 


oc = 10° 


(X = 20 ° 




Mach Contours 


Error Controlled Aero Database 




(X = 0 ° 


oc = 10° 


(X = 20 ° 



Mach Contours 




Error Controlled Aero Database 




(X = 0 ° 


oc = 10° 


(X = 20 ° 







Error Controlled Aero Database 




(X = 0 ° 


oc = 10° 


(X = 20 ° 



Mach Contours 



Error Controlled Aero Database 



(X = 0 ° 


oc = 10° 



Mach Contours 


Error Controlled Aero Database 



Moo = 1.6 

(X = 0 ° 


oc = 10° 


(X = 20 ° 



Mach Contours 



Error Controlled Aero Database 




(X = 0 ° 


oc = 10° 


(X = 20 ° 





Mach Contours 



Error Controlled Aero Database 



Normal Force 
Coefficient 
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Error Controlled Aero Database 



Axial Force 
Coefficient 



Error Controlled Aero Database 


Pitch moment 
about nose 0 
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Error Controlled Aero Database 




Error Controlled Aero Database 
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Error Controlled Aero Database 




Error Controlled Aero Database 
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• • Functional 

Corrected Functional 
G O Error Bound Estimate 
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Launch Abort System 
Tower Jettison Database 



Mach -0.75 
6 sec 

-5300 ft altitude, -3300 ft downrange 


Mach -0.75 



Mach -0.3 
15 sec 



Mach -0.2 
20 sec 


ACM trims LAV at -0 deg alpha during 
ACM coast phase. Slows down from 
Mach 0.8 to 0.2 in -10 seconds and 
drops in qbarfrom 800 to 100 psf 


ACM holds LAV at -25 deg 
alpha during AM burnout. 


4 


Mach -0.3 
1 sec 


ACM + Canards turn LAV to 
heat shield forward 



ACM initially pushes 
LAV to -25 deg alpha 
during AM burn 


8 Abort Control 


4 Abort 


Mach 0.2 
-21.1 sec 


LAS is jettisoned with ACM 
still burning. CM in free flight 
for -3 seconds until drogue 
chute deployment. 



Motors (AC Ms) Motors (AMs) 


Mach 0 
0 sec 

Ignition 






4 Separation 
Motors (SMs) 


C.G. 


Launch Abort System Tower Jettison 


• Objective: Analyze aerodynamic forces and moments during LAS tower 
jettison 


• Include effects of jettison motor firing with translation and rotation of the Crew 
Module (CM) relative to the Launch Abort Module 



► N 



Geometry Configuration 


Four configuration parameters for 
CM position and orientation: 

Ax, Ay, Az, 0 


A ZCM 


A XCM 


Database Cases 




• Run Conditions: 

Moo = {0.5, 0.7, 0.9, 1.1, 1.3, 1.6} 

« = {155°, 160°, 165°, 170°, 175°, 180°} 
P = {0°, 5°} 

CMAx, CMAy, CMAz, CMA0 


Jettison Motor Plume 
Conditions: 

• JM on and JM off 

• Scale thrust for altitude 


Trajectory: maintain ~ Const, 


1200 


.* 




A XCM 


Functional and Adaptation Strategy 


• Challenging simulations 

► Complex, detailed geometry 

► Bodies in close proximity 

► Strong upstream-firing jets 

► Shocks and wakes 



• Output functional: J = ( 0 . 8 | Gv | + 02\Ca\)cm + (| Ov | + 0A\Ca\)lam 


• Solution technique is a compromise since most of these cases are unsteady, 
and need high resolution 

► Want the best answer as cheaply as possible 

• Adaptation follows “worst-things-first” strategy (AIAA 2008-0725) 

► Refine largest contributors to output error first 


Capture forces on CM 


Capture forces on tower 


Initial Mesh 



• Background mesh essentially symmetric and coarse to avoid biasing solution 


p_> 


3000 Cells 




Example Final Mesh (1 0 Adapt Cycles) 


• Output functional: J = ( 0 . 8 | Gv | + 02\Ca\)cm + (| CV | + 0A\Ca\)lam 

• Mao = 1 .1 , « = 1 60° with the CM @ (Ax, Ay, Az, 1 0°) 



Example Final Solution (1 0 Adapt Cycles) 


• Output functional: J = (0.8|Gv| + 02\Ca\)cm + (| CV | + 0A\Ca\)lam 


• Moo = 1 .1 , « = 1 60° with the CM @ (Ax, Ay, Az, 1 0°) 



Mach 

ontours 


2.95 M 
Cells 



Convergence of Aerodynamic Coefficients 


• Output functional: J = ( 0 . 8 | C / v | + 0.2\Ca\)cm + (| CV | + QA\Ca\)lam 


• M oo = 1 .1 , a = 1 60° with the CM @ (Ax, Ay, Az, 1 0°) 


• Convergence of 
forces and moms, 
on CM with mesh 
refinement 


t 1 — f — i 1 j 1 1 f 1 j| 1 1 1 i 1 ii 1 j 1 1 r 
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0 500 1000 1500 2000 

Multigrid Cycles 


Iteration 

averages 


Comparison with Unsteady Simulation 


• How do these iteration averages compare with unsteady simulation? 

• Performed comparisons at Mach 0.7, & 1.1 

• “Best unsteady mesh” generated by 1 additional refinement of steady mesh, 
using low threshold to “fill out” adaptation regions 

• Example case: = 1 .1 , a = 1 60° with the CM @ (Ax, Ay, Az, 1 0°) 


Steady mesh & Cp, -2.95 M Cells Unsteady mesh & Cp, -5.2 M Cells 




Comparison with Unsteady Simulation 



• How do these iteration averages compare with unsteady simulation? 

• Performed comparisons at Mach 0.7, & 1.1 

• “Best unsteady mesh” generated by 1 additional refinement of steady mesh, 
using low threshold to “fill out” adaptation regions 

• Example case: = 1 .1 , a = 1 60° with the CM @ (Ax, Ay, Az, 1 0°) 



Steady mesh iso-Mach, -2.95 M Cells 



Unsteady mesh iso-Mach, -5.2 M Cells 



Model Axis Force and Moment Coef. 


Comparison with Unsteady Simulation 


M oo = 1 .1 , « = 1 60° with the CM @ (Ax, Ay, Az, 1 0°) 

Steady, -2.95 M Cells Unsteady, -5.2 M Cells 



• Agreement to third significant digit 

• Differences in averages are same size as differences due to averaging window 

• Similar results for other components at both Mach 0.7 & Mach 1 .1 


Plume Shape 


To iso-surface showing approximate plume shape 



Moo = 1.3, a = 160° 


Adaptation chases only those 
regions of plumes that effect 
loads 


68 
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Comparison of JM On and Off Flowfields 




> M cells 


= 1 .1 , a = 1 55°, 2.85 M cells 
Surface C P 


Mach Contours 


JM plumes move bow shock ~1 8 ft upstream 


JM Power On, LAM Only 


Database Samples 


Loads on CM and LAM in proximity 

(Ax, 0.0, 0.0, 0°) 


Loads on CM 
Loads on LAM 





Cm (Pitch) C|\lormal Cm (Pitch) ClMormal 





Database Samples 


Loads on CM and LAM in proximity 

(Ax, 0.0, 0.0, 0°) 


Loads on CM 
Loads on LAM 





Ci\lormal 





Database Samples 


Loads on CM and LAM in proximity 

(Ax, 0.0, 0.0, 0°) 


Loads on CM 
Loads on LAM 




Cl\lormal 





Database Samples 


Loads on CM and LAM in proximity 

(Ax, 0.0, 0.0, 0°) 


Loads on CM 
Loads on LAM 







Database Samples 


Loads on CM and LAM in proximity 

(Ax, 0.0, -Az, 0°) 


Loads on CM 
Loads on LAM 





Cm (Pitch) C|\l ormal Cm (Pitch) CNormal 






Database Samples 


Loads on CM and LAM in proximity 

(Ax, 0.0, Az, 0°) 


Loads on CM 
Loads on LAM 





Cm (Pitch) C|\l ormal Cm (Pitch) CNormal 






Summary 

Presented a reliable and efficient approach for error estimates 
and mesh refinement of complex geometry problems 

1 . Handles complex geometry problems in an automatic fashion 

2. Tolerant of coarse initial meshes 

3. Behavior of functional, correction, and error estimate provide an 
indication of errors due to lack-of-convergence in steady simulations 

It is our best mesh generator ... refinement 
^ complements and surpasses expert 
knowledge 

=|> Allows users to focus on data validation and 
analysis instead of mesh generation 


Present and Future Work 


Sonic-boom applications (Mathias Winzter, AIAA 2008-6593) 
Address unsteadiness issues in difficult cases 

- Affordable mesh refinement and error bound for “mildly” 
unsteady flow 

- Formal unsteady adjoint development 

Control accuracy of objective functions in optimization studies 
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